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ABSTRACT 

Context. Double-lined, detached eclipsing binaries are our main source for accurate stellar masses and radii. In this 
paper we focus on the 1.15-1.70 Mq interval where convective core overshoot is gradually ramped up in theoretical 
evolutionary models. 

Aims. We aim to determine absolute dimensions and abundances for the F-type detached eclipsing binary BK Peg, and 
to perform a detailed comparison with results from recent stellar evolutionary models, including a sample of previously 
studied systems with accurate parameters. 

Methods, uvby light curves and uvby/3 standard photometry were obtained with the Stromgren Automatic Telescope, 
ESO, La Silla, and high-resolution spectra were acquired with the FIES spectrograph at the Nordic Optical Telescope, 
La Palma. 

Results. The 5?49 period orbit of BKPeg is slightly eccentric (e = 0.053). The two components are quite different with 
masses and radii of (1.414±0.007 Mq, 1.988± 0.008 Rq) and (1.257±0.005 Mq, 1.474±0.017 Rq), respectively. The 
measured rotational velocities are 16.6±0.2 (primary) and 13.4±0.2 (secondary) kms -1 . For the secondary component 
this corresponds to (pseudo) synchronous rotation, whereas the primary component seems to rotate at a slightly lower 
rate. We derive an iron abundance of [Fe/H] = —0.12 ± 0.07 and similar abundances for Si, Ca, Sc, Ti, Cr and Ni. The 
stars have evolved to the upper half of the main-sequence band. Yonsei-Yale and Victoria-Regina evolutionary models 
for the observed metal abundance reproduce BK Peg at ages of 2.75 and 2.50 Gyr, respectively, but tend to predict 
a lower age for the more massive primary component than for the secondary. We find the same age trend for three 
other upper main-sequence systems in a sample of well studied eclipsing binaries with components in the 1.15-1.70 Mq 
range. We also find that the Yonsei-Yale models systematically predict higher ages than the Victoria-Regina models. 
The sample includes BW Aqr, and as a supplement we have determined a [Fe/H] abundance of —0.07 ± 0.11 for this 
late F-type binary. 

Conclusions. We propose to use BKPeg, BWAqr, and other well-studied 1.15-1.70 Mq eclipsing binaries to fine-tune 
convective core overshoot, diffusion, and possibly other ingredients of modern theoretical evolutionary models. 

Key words. Stars: evolution - Stars: fundamental parameters - Stars: binaries: eclipsing - Stars: individual: BKPeg, 
BWAqr - Techniques: photometric - Techniques: spectroscopic 



^ ' 1. Introduction al. (|2008bl hereafter CTB08). One of the troublesome ingre- 
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to 1% or better (Torres et al. I2009|) . and they also pro- 
vide stringent tests of various aspects of stellar evolution- xhe literature on the existence and calibration of core 

ary models. For this purpose, well-established abundance overshoot is extensive, and here we only draw attention to a 
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Fig. 1. y light curve and b — y and u — b colour curves (instrumental system) for BK Peg. 



0.18 (in units of the pressure scale height). The question of 
mass-dependence of the degree of core overshoot has - again 
based on dEB samples - been addressed by e.g. Ribas et al. 
([20"0"0"f and Claret pOOTj) . but they arrive at very different 
conclusions. 

Most recent grids of stellar evolutionary calculations in- 
clude core overshoot, but the recipes in terms of mass and 
abundance dependence are somewhat different. We refer to 
Sect. GO for details. 

Our main motivation to undertake a study of BK Peg 
has been that a) it has evolved to the upper part of the 
main-sequence band and is therefore well-suited for core 
overshoot tests, but published dimensions are not of suffi- 
cient quality, and b) few similar well-studied systems are 
known. Below we present absolute dimensions and abun- 
dances based on new uvby light curves and high-resolution 
spectra and compare BK Peg and other similar systems 
with Yonsei-Yale and Victoria-Regina stellar evolutionary 
models. 

In Appendix [X] we present a spectroscopic abundance 
analysis for BWAqr as supplement to the study of this 
binary by Clausen (|199ip . 

2. BK Peg 

BKPeg (BD+25 5003, m v = 9.98, Sp. type F8, P = 
5^49), is a well detached, double- lined eclipsing binary with 
1.41 and 1.26 Mq components in a slightly eccentric (e = 



0.0053) orbit. It is unusual in the sense that the more mas- 
sive, larger, and more luminous component is slightly cooler 
than the other component. This is due to evolution, where 
the more massive component has evolved to the upper part 
of the main-sequence band. 

The eclipsing nature of BKPeg was discovered by 
Hoffmeister ([iggTji . and Lause PM51 HMTf visual ly ob- 
served several times of minima. Popper & Dumont (|1977[) 
obtained BV light curves at the Palomar and Kitt Peak ob- 
servatories, which were later analysed by Popper & Etzcl 
(1981). Preliminary absolute dimensions were included by 
Popper in his critical review of stellar masses (Popper 
1980), and soon after he presented spectroscopic elements 
and improved absolute dimensions (Popper 1983). He 
reached component masses accurate to about 1% and radii 
accurate to 4% (primary) and 7% (secondary). Demircan et 
al. (|1994|) presented photoelectric U BV light curves and an 
improved ephemeris, as well as absolute dimensions, which, 
within the uncertainties, agree with those by Popper; see 
also Popper & Etzel (|1995|) for clarification on terminology. 
We refer to the more massive, larger but cooler component 
as the primary (p) component, which, for the ephemeris we 
adopt (Eq. Q}, is eclipsed at phase 0.0. 

3. Photometry 

Below, we present the new photometric material for BK Peg 
and refer to Clausen et al. (|20011 hereafter CHO01) for 



J.V. Clausen et aL: BKPeg 
Table 1. Photometric data for BKPeg and the comparison stars. 
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a Popper (Tl983"j) . b Harlan & Taylor (T1970)) . c Heard (Tl956j) . d Harlan (Tl969)l 

NOTE 1: References are: C10 = This paper. J96 = Jordi et al. (Tl996)> . 083 = Olsen (fl983"j) . 094 = Olsen |Tl994)l . P69 = Perry 
GSSHJ). 

NOTE 2: For BKPeg, the uvby/3 information by C10 is the mean value outside eclipses (first line) and during the central part of 
secondary eclipse (second line), where 98% (y) of the light of the secondary component is eclipsed. 

NOTE 3: N is the total number of observations used to form the mean values, and a is the rms error (per observation) in mmag. 



further details on observation and reduction procedures, 
and determination of times of minima. 

3.1. Light curves for BK Peg 

The differential uvby light curves of BK Peg were observed 
at the Stromgren Automatic Telescope (SAT) at ESO, La 
Silla with its 6-channel uvbyft photometer on 66 nights 
between October 2000 and September 2003 (JD2451828- 
2452910). They contain 384 points per band with most 
phases covered at least twice. The observations were done 
through an 18 arcsec diameter circular diaphragm at air- 
masses between 1.8 and 2.2. HD 222249, HD 222391, and 
HD 223323 - all within a few degrees of BK Peg on the sky - 
were used as comparison stars and were all found to be con- 
stant within a few mmag; see Table [1] The light curves are 
calculated relative to HD 223323, but all comparison star 
observations were used, shifting them first to the same light 
level. HD 223323 has later been found to be a double-lined 
spectroscopic binary with an estimated orbital inclination 
of about 63°; the orbital period is 1175 days and the orbital 
eccentricity is 0.6 (Griffin 2007). The average accuracy per 
light curve point is about 5 mmag (ybv) and 8 mmag (it). 
The light curves (Table 13) will only be available in elec- 
tronic form. 

As seen from Fig.[T] BK Peg is well detached with nearly 
identical eclipse depths of about 0.5 mag. In the y, b, and v 
bands, the primary eclipse at phase 0.0, which is a transit, 
is slightly deeper than the secondary eclipse (almost total 
occultation) , which occurs at phase 0.4976. In the u band 
the secondary eclipse is, however, the slightly deeper one. 

3.2. Standard photometry for BK Peg 

Standard uvby/3 indices for BK Peg and the three compari- 
son stars, observed and derived as described by CHO01, are 
presented in Table[T] As seen, the indices are based on many 



observations and their precision is high. For comparison, we 
have included published photometry from other sources. In 
general, the agreement is good, but in some cases differ- 
ences are larger than the quoted errors; we have used the 
new results for the analysis of BK Peg. 

3.3. Times of minima and ephemeris for BK Peg 

Three times of secondary minimum, but none of primary, 
have been determined from the uvby light curve observa- 
tions. They are listed in Table [5] together with available 
measured times. From separate weighted least squares fit to 
the times of primary and secondary minima, respectively, 
we derive the linear ephemeris given in Eq. [1] 

Min 1 = 2450706.46968 + 5 d 48991046 x E m 
±38 ±56 1 ' 

Within errors, the two types of minima yield identical peri- 
ods, and the new ephemeris is in good agreement with that 
by Demircan et al. (|1994p ; see also Kreiner et al. (|200ip and 
Kreiner ([StM)! 1 ! . 

4. Spectroscopy 

In order to perform abundance determinations and also 
improve the spectroscopic elements by Popper (1983), 
we have obtained 13 high-resolution (R = 45000) spec- 
tra with the FIES fibre echelle spectrograph at Nordic 
Optical Telescope, La Palma during five consecutive nights 
in August 2007; see Table [3] For the basic reduction of the 
spectra, we have applied the IRAF based FIEStool pack- 
aged Subsequently, dedicated IDI0 programs were applied 

|http : // www . as . ap . krakow . pi/ ephem 

2 see http://www.not.iac.es for details on FIES and 
FI EStool. 

3 http : //www. ittvis . com/idl/index . asp 



4 J.V. Clausen et al.: BKPeg 

Table 2. Times of primary (P) and secondary (S) minima for BKPeg. 
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a O-C values (P) and phases (S) are calculated for the ephemeris given in Eq. [T] 
b Observing methods are: V = visual; PE = photoelectric; CCD = CCD. 

c References are: A03 = Ak et al. (2003}. B87 = Braune & Hiibscher ([T987)l . C10 - This paper. D 77 = D worak [Tl977)l . D94 = 
Demircan et al. ([T5M)> . DP I = D iethelm ipTJOTjl . H07 = Hiibscher & Walter (JUT}. L35 = Lause 1 (1535)1 . L37 = Lause 1(1537)1 . 
P83 = Popper & Dumont (I1977p . redetermined. 



to remove cosmic ray events and other defects, and for nor- 
malisation of the individual orders. For each order, only the 
central part with acceptable signal-to-noise ratios was kept 
for further analysis. 

The radial velocities for BK Peg were measured from 40 
useful orders of the 13 FIES spectra. We applied the broad- 
ening function (BF) formalism (Rucinski 115551 [20021 I2TJ04")) . 
using v sin i = synthetic templates matching the ef- 
fective temperature, log(g), and metal abundance of the 
components of BKPeg. They were calculated with the 
bssynth tool, which applies the SYNTH software (Valenti 
& PiskunovdHSl) and modified ATLAS9 models (Heiter et 
al. [5002). Since the components have nearly identical tem- 
peratures, the two templates are very similar and lead to 
practically identical results. As described by e.g. Kaluzny 
et al. (|2006p . the projected rotational velocities vsmi of 
the components and (monochromatic) light /luminosity ra- 
tios between them can also be obtained from analyses of 
the BFs. 



For each spectrum, BFs were calculated for each of the 
selected orders, and a mean BF was then calculated to- 
gether with weights for each order based on the root mean 
square deviation of the individual BFs from the mean BF. 
The final BF is the weighted average for the selected or- 
ders. The radial velocites, the vsini's, and the light ratio 
were derived by fitting a rotational profile for both stellar 
components, convolved with a Doppler profile correspond- 
ing to the instrumental resolution, to the final BF for each 
observed spectrum. 

The radial velocities are listed in Table 2) The final val- 
ues of v sin i and light ratio were calculated as the mean 
values for the 13 spectra, with errors estimated from the de- 
viations from spectrum to spectrum. For the primary and 
secondary components of BK Peg, we obtain mean rota- 
tional velocities of 16.6±0.2 and 13.4±0.2 kms -1 , respec- 
tively. For the light ratio we find L s /L p = 0.57 ± 0.02. 

In addition, we have determined the light ratio between 
the components by directly comparing the FIES spectra 
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Table 3. Log of the FIES observations of BKPeg. 



HJD-2 400 000^ phase t exp ^ S/N£ 

60 
55 
20 
70 
85 
85 
100 
55 
60 
75 
80 
65 
75 



a Refers to mid-exposure 
b Exposure time in seconds 

c Signal-to-noise ratio measured around 6070 A 

Table 4. Radial velocities of BK Peg and residuals from 
the final spectroscopic orbit presented in Table [5J 



HJD Phase RV P RV S (O - C) p {0-C)„ 

-2 400 000 kmr 1 kms -1 kms -1 kms -1 

0.28 
0.38 
0.09 
0.38 
0.33 
0.42 
0.10 
0.44 
0.47 
0.45 
0.50 
0.59 
0.80 



and synthetic binary spectra, calculated for a range of 
luminosity ratios between the components. Adopting the 
temperatures, surface gravities, rotational velocities, and 
metallicities listed in Table ^Ul and using several spectral 
orders covering 5300-5800 A, we obtain the best line fits 
for a light ratio of 0.57 ± 0.04. As expected, since the com- 
ponents have nearly identical temperatures, we find no sig- 
nificant wavelength dependence of the spectroscopic light 
ratio, even if a broader wavelength region is used. 

5. Photometric elements 

Since BKPeg is well-detached, the photometric ele- 
ments have been determined from jktebofQ analyses 
(Southworth et al. 2004a 2004b) of the uvby light curves. 
The underlying simple Nelson-Davis-Etzel binary model 
(Nelson fc Davi s [W2l Etzel [TMT1 Popper & Etzel fTMTl 
Martynov I1973[) represents the deformed stars as biaxial 
ellipsoids and applies a simple bolometric reflection model. 
We refer to CTB08 for details on the general approach ap- 
plied. In tables and text, we use the following symbols: i 
orbital inclination; e eccentricity of orbit; lu longitude of 

4 http : //www. astro .keele . ac.uk/~jkt/ 



periastron; r relative radius; k = r s /r p ; u linear limb dark- 
ening coefficient; y gravity darkening coefficient; J central 
surface brightness; L luminosity; T c g effective temperature. 

The mass ratio between the components was kept at 
the spectroscopic value, see Sect. [6l The simple built-in 
bolometric reflection model was used, linear limb darken- 
ing coefficients by Van Hamme (|1993p and Claret (2000) 
were applied, and gravity darkening coefficients correspond- 
ing to radiative atmospheres were adopted. Identical coef- 
ficients were used for the two components, since their ef- 
fective temperatures and surface gravities are sufficiently 
identical. Effective temperatures determined from the stan- 
dard uvby and JHK S indices outside eclipses are listed in 
Table [5) As seen, the results from the different calibrations 
agree well; we have adopted the temperature based on the 
Holmberg et al. (2007) calibration. 

Solutions for BK Peg, based on Van Hamme limb dark- 
ening coefficients, are presented in Table [51 and O—C resid- 
uals of the y observations from the theoretical light curve 
are shown in Fig. [21 As seen, the results from the four bands 
agree well. Changing to Claret (|2000p limb darkening coef- 
ficients, which are 0.07-0.09 higher, increases the radius of 
the primary component by only 0.4%, whereas that of the 
secondary component is increased by 1.5%. This is linked 
to a 1% larger k and a 29% smaller esin(w), reducing e by 
10%. Limb darkening coefficients determined from the light 
curves reproduce those by Van Hamme better than those 
by Claret, but have uncertainties of about ±0.12. Including 
non-linear limb darkening (logarithmic or square-root law) 
has no significant effect on the photometric elements. 

The adopted photometric elements listed in Table [7] are 
the weighted mean values of the jktebop solutions adopt- 
ing the linear limb darkening coefficients by Van Hamme. 
Realistic errors, based on 10 000 Monte Carlo simulations in 
each band and on comparison between the uvby solutions, 
have been assigned. The Monte Carlo simulations include 
random variations within ±0.07 of the linear limb darken- 
ing coefficients. As seen, r p becomes more accurate than r s . 
This is because it correlates less with k, probably related 
to the secondary eclipse being nearly total for the adopted 
elements. It should be noted that the ybv luminosity ra- 
tios from the light curve solutions agree very well with the 
spectroscopic light ratio (Sect. 01 . 

For comparison, Popper & Etzel (|1981[) obtained r p = 
0.108 ± 0.04 and r s = 0.086 ± 0.006, assuming esin(u;) = 
and adopting k — 0.80 ± 0.05 in order to reproduce a mean 
ratio of 0.69 ±0.03 between selected secondary and primary 
lines, measured on photographic spectra. This ratio is, how- 
ever, much higher than the light ratio we derive from the 
FIES spectra, leading to a higher k. Popper & Etzel also 
obtain a somewhat larger value for r p + r s . The relative 
radii presented by Demircan et al. (1994) are close to those 
by Popper & Etzel. 

In conclusion, the new photometric elements derived 
from analyses of the uvby light curves are significantly more 
accurate than previous determinations. We find that the 
secondary eclipse is almost total, with 98% of the y light 
of the secondary component eclipsed, whereas about 57% 
of the y light from the primary component is eclipsed at 
phase 0.0. 
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54333.45954 


0.6647 


60.398 


-83.802 


0.25 


54333.46395 


0.6655 


60.694 


-83.936 


0.35 


54334.42028 


0.8397 


58.578 


-81.562 


-0.50 


54335.52799 


0.0414 


-28.144 


15.916 


-0.19 


54335.53872 


0.0434 


-28.875 


16.995 


0.02 


54335.62203 


0.0586 


-36.037 


24.823 


0.05 


54335.72881 


0.0780 


-44.698 


34.832 


0.22 


54336.41869 


0.2037 


-83.237 


77.803 


0.15 


54336.43503 


0.2067 


-83.671 


78.279 


0.11 


54336.52436 


0.2229 


-85.447 


80.223 


0.03 


54336.74168 


0.2625 


-86.151 


81.269 


-0.01 


54337.44772 


0.3911 


-56.696 


48.134 


-0.20 


54337.58944 


0.4169 


-46.161 


36.119 


-0.28 
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Table 5. Effective temperatures (K) for the combined light of BKPeg. 



A v 


[Fe/H] 


(b - y) c (V - J) 


(V - H) {V - K s )o A96 


H07 


RM05 


M06 


0.188 


-0.12 


0.405 0.320 1.208 


1.458 1.564 6280 


6270 


6270/6255/6130/6265 


6290 



Note 1: Av is the adopted visual interstellar absorption. The Vo magnitude and the (b — y)o and Co indices are based on the 
out-of-eclipse uvby standard indices from Table[T] The 2MASS observations (J,H,K a ) were obtained at phase 0.879. 
Note 2: Referen ces ar e: A96 = Alonso et al. Q996J. H07 = Holmberg et al. (12007)) . RM05 = Ramirez & Melendez (I20051 . M06 
= Masana et al. ipOTTp]) . 

Note 3: The results from A96 are based on their uvby calibration, those from RM05 on their uvby, (V — J), (V — H), and (V — K s ) 
calibrations (in that order); the calibration by M06 is for (V — K s ). 



Table 6. Photometric solutions for BK Peg from the JK- 
tebop code. 





V 


b 


V 


u 


i n 


88.02 


87.91 


87.99 


87.90 


±5 


±5 


±4 


±7 


ecosu; 


-0.00364 


-0.00376 


-0.00363 


-0.00374 




±6 


±7 


±6 


±10 


esin uj 


0.00283 


0.00563 


0.00438 


0.00059 




±181 


±207 


±183 


±312 




0.0046 


0.0068 


0.0057 


0.0038 


I n 


142.1 


123.8 


129.6 


170.0 


r p 


0.1092 


0.1096 


0.1092 


0.1100 


r a 


0.0806 


0.0819 


0.0807 


0.0813 


k 


0.7379 


0.7474 


0.7394 


0.7391 




±50 


±68 


±46 


±97 


r p + r 3 


0.1898 


0.1915 


0.1900 


0.1913 




±5 


±5 


±4 


±7 


Up — Us 


0.55 


0.64 


0.72 


0.70 


Vp = Vs 


1.09 


1.24 


1.42 


1.68 


J s 1 Jp 


1.0444 


1.0447 


1.0569 


1.1016 


±31 


±37 


±34 


±60 


L s / L p 


0.5670 


0.5817 


0.5756 


0.5990 


a (mmag.) 


4.8 


5.2 


4.7 


7.9 



0.02 
'So 0.01 

■3 0.00 

u 

i -0.01 
-0.02 



~~ Primary eclipse 




+ + 1 +' +++ + + ++ + 





-0.03 -0.02 -0.01 



0.00 
Phase 



0.01 0.02 



0.02 
0.01 



O -0.01 
-0.02 



Secondary eclipse 

+ 

++++ + 


+ + + 




++ + + ++ + 

+ + + 



0.47 0.48 0.49 0.50 0.51 0.52 

Phase 



0.02 
01 



O -0.01 
-0.02 



Note 1: Linear limb darkening coefficients by Van Hammc 
( 1993!) were adopted, a mass ratio of 0.89 was assumed, and 
phase shift and mag. normalization were included as free pa- 
rameters 

Note 2: The errors quoted for the free parameters are the 
formal errors determined from the iterative least squares so- 
lution procedure 



0.02 

^ 001 
— 0.00 
° -0.01 

-0.02 



Table 7. Adopted photometric elements for BK Peg. 



i 

e 

r p 
r s 
r p 
k 



+ r 3 



87? 96 ±0? 14 
0.0053 ±0.0013 
138? 7 ± 19? 
0.1094 ±0.0004 
0.0811 ±0.0009 
0.1906 ±0.0012 
0.741 ± 0.008 



J s I Jp 



y 

1.044 
±4 
0.572 
±11 



b 

1.048 
±4 
0.574 
±14 



1.057 
±4 
0.578 
±10 



1.095 

±6 
0.599 
±15 



Note: The individual flux and luminosity ratios are based on 
the mean stellar and orbital parameters 



Fig. 2. (O—C) residuals of the BK Peg y-band observations 
from the theoretical light curve computed for the photomet- 
ric elements given in Table [5] 



6. Spectroscopic elements 

Spectroscopic orbits have been derived from a re-analysis 
of the radial velocities by Popper (1983) and an analysis of 
the new radial velocities listed in Table |U We have used the 
method of Lehman-Filhes implemented in the SBOF0 pro- 
gram (Etzel[2004), which is a modified and expanded ver- 
sion of an earlier code by Wolfe, Horak & Storer (1967J). The 
orbital period P was fixed at the ephemeris value (Eq. [1]), 
and the eccentricity e and longitude of periastron uj to the 
results from the photometric analysis (Table [7]) . The radial 



5 Spectroscopic Binary Orbit Program, 
http://mintaka.sdsu.edu/faculty/etzel/ 
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Table 8. Spectroscopic orbital solutions for BKPeg. 



Parameter: 


Popper 


FIES 




(re- analysis) 


(adopted) 



Adjusted quantities: 
K p (kms^ 1 ) 
K s (kms" 1 ) 
7 P (kms -1 ) 
7 S (kms -1 ) 

Adopted quantities: 
P (days) 

T (HJD-2 400 000) a 



u (°) 

Derived quantities: 
M p sin 3 i (Mq) 
M s sin 3 i (Mq) 
q = M s /M p 
asini (Rq) 

Other quantities 

pertaining to the fit: 

N obs (p/ s) 

Time span (days) 

a p h (kms -1 ) 

a s h (kms" 1 ) 



79.14 ±0.33 
88.86 ± 0.56 
-8.51 ±0.30 
-8.35 ±0.51 

5.48991046 
50706.46968 
0.0053 
138.7 

1.427 ±0.019 
1.271 ±0.013 
0.891 ±0.007 
18.231 ±0.071 



24/23 
1374 
1.36 
2.30 



78.77 ±0.11 
88.59 ±0.21 
-7.39 ±0.08 
-7.20 ±0.15 

5.48991046 
50706.46968 
0.0053 
138.7 

1.411 ±0.007 
1.255 ±0.005 
0.889 ± 0.002 
18.161 ±0.026 



13/13 
4 

0.26 
0.50 



a Time of central primary eclipse 

b Standard deviation of a single radial velocity 



velocities of the components were analysed independently 
(SB1 solutions). 

The spectroscopic elements are presented in Table |5] 
The semiamplitudes (K p ,K s ), and their uncertainties, ob- 
tained from Popper's velocities are identical to his results, 
even though he assumed the orbit to be circular and used 
an older ephemeris. 

As seen, significantly more accurate semiamplitudes are 
derived from the FIES velocities. They are slightly smaller 
than those from Popper's velocities, but within errors the 
results agree. Including e and/or u) as free parameters for- 
mally improves the solution but does not alter the semi- 
amplitudes. The number of velocities is, however, too small 
for reliable spectroscopic determination of e and lu. The 
double-lined (SB2) solutions agree perfectly with the single- 
lined solutions listed in Table [8] Also, spectroscopic el- 
ements determined as part of the spectral disentangling 
(Sect. [7]) are identical. We notice that the new system veloc- 
ities (7 P ,7s) differ by about 1 kms -1 from Popper's results. 
This is probably due to radial velocity zero point differ- 
ences. Our velocities are tied to the ThAr exposures taken 
before and/or after each target exposure. Standard star ob- 
servations normally agree to within 0.1-0.2 kms -1 . 



7. Chemical abundances 

For abundance analyses, we have disentangled the FIES 
spectra of BK Peg in order to extract the individual com- 
ponent spectra. We have applied the disentangling method 
introduced by Simon & Sturm ( 1994J) and a revised version 
of the corresponding original code developed by E. Sturm. 
It assumes a constant light level, but since BK Peg is con- 
stant to within 0.5% outside of eclipses this is of no concern. 
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Fig. 

line: 



3. Spectroscopic orbital solution for BKPeg (solid 
primary; dashed line: secondary) and radial velocities 



(filled circles: primary; open circles: secondary). The dot- 
ted line (upper panel) represents the center-of-mass veloc- 
ity of the system. Phase 0.0 corresponds to central primary 
eclipse. 

Twenty- two orders, covering 5160-6450 A (with a few gaps) 
were selected and disentangled individually. The orbital el- 
ements were fixed at the adopted values (Table [SJ and very 
slightly wavelength dependent light ratios matching the re- 
sults from the light curve analyses (Table [7]) were adopted. 
Around 6070 A, the signal-to-noise ratios of the resulting 
component spectra are 160 (primary) and 80 (secondary). 
A 40 A region, centred at 6070 A, is shown in Fig. [4] 

The basic approach followed in the abundance analyses 
is described by CTB08. We used the versatile VWA tool, 
which applies the SYNTH software (Valenti & Piskunov 
1996) to compute synthetic spectra. We refer to Bruntt et 
al. (PUU41I2"0TJ81) and Bruntt (|3D05|l for a detailed description 
of VWA. Atmosphere models were interpolated from the 
recent grid of MARCS model atmospheres (Gustafsson et 
al. I2008p . which adopt the solar composition by Grevesse 
et al. (|2007p . Atomic line data are from the Vienna Atomic 
Line Database (VALD; Kupka et al. I1999p . but in order to 
derive abundances relative to the Sun, log(gf) values have 
been adjusted in such a way that each measured line in the 
Wallace et al. (1998) Solar atlas reproduces the atmospheric 
abundances by Grevesse et al. (20071). 

The abundance results derived from all useful lines with 
equivalent widths between 10 and 100 mA are presented 
in Table [5] The equivalent widths measured in the dis- 
entangled spectra are listed in Tables 14 (primary) and 
15 (secondary), which will only be available in electronic 
form. The surface gravities and observed rotational veloci- 
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Fig. 4. A 40 A region centred at 6070 A of the disentangled spectra of the components of BK Peg. Lines identified by a 
red line were not used for the abundance analysis. 



Table 9. Abundances ([El./H]) for the primary and sec- 
ondary components of BK Peg. 





P 


rimary 




Secondary 




Ion 


[El./H] 


rms 


N a 


[El./H] 


rms 


N a 


Si I 


-0.19 


0.05 


12 


-0.12 


0.09 


10 


Cal 


-0.05 


0.05 


4 


0.00 


0.07 


6 


Sell 


-0.14 


0.09 


3 


-0.18 


0.09 


3 


Til 








-0.17 


0.06 


3 


Tin 


-0.07 


0.09 


3 


0.01 


0.18 


3 


Cri 


-0.20 


0.03 


3 


-0.19 


0.03 


3 


Cm 


-0.23 


0.05 


3 


-0.04 


0.03 


4 


Fei 


-0.12 


0.07 


78 


-0.12 


0.09 


70 


Fen 


-0.15 


0.09 


10 


-0.11 


0.11 


12 


Nil 


-0.23 


0.08 


6 


-0.24 


0.13 


11 



Number of lines used per ion 



ties listed in Table [10] were adopted, whereas the effective 
temperatures and microturbulence velocities were tuned 
until Fe I abundances were independent of line equivalent 
widths and excitation potentials. The resulting tempera- 
tures are 6365 ± 50 K (primary) and 6385 ± 50 K (sec- 
ondary). From a study of 10 F-K type stars with interfero- 
metrically and spectroscopically determined effective tem- 
peratures, Bruntt et al. (|2010p find a systematic offset of 
40 K, which should be subtracted. The corrected spectro- 
scopic temperatures are still slightly higher than derived 
from the uvby indices ( Table [TU]) but agree within errors. 

Microturbulence velocities of 1.55 ± 0.25 kms" 1 (pri- 
mary) and 1.22 ± 0.25 kms -1 (secondary) were obtained. 
The calibration by Edvardsson et al. () 1993)) predicts higher 
values of 2.20±0.31 kms' 1 (primary) and 1.95±0.31 kms -1 
(secondary); the difference in microturbulence will be dis- 
cussed by Bruntt et al. (|2010p . 

As seen, a robust [Fe/H] is obtained, with nearly identi- 
cal results from Fe I and Fe n lines of both components. 



The mean value from all measured Fe lines is [Fe/H] = 
—0.12 ± 0.01 (rms of mean). Changing the model tempera- 
tures by ±50 K modifies [Fe/H] from the Fe I lines by about 
±0.05 dex, whereas almost no effect is seen for Fe II lines. If 
0.25 kms -1 higher microturbulence velocities are adopted, 
[Fe/H] decreases by about 0.04 dex for both neutral and ion- 
ized lines. Taking these contributions to the uncertainties 
into account, we adopt [Fe/H] = -0.12 ± 0.07 for BKPeg. 
In general, we find similar relative abundances for the other 
ions listed in Table [9j including the a-elements Si, Ca, and 
Ti. 

As an addition to the spectroscopic abundance anal- 
ysis, we have also calculated metal abundances from the 
de-reddened uvby indices for the individual components 
(Table [TU]) and the calibration by Holmberg et al. (|2007p . 
The results are: [Fe/H] = -0.07 ± 0.18 (primary) and 
[Fe/H] = —0.14 ± 0.20 (secondary). Within errors they 
agree with those from the spectroscopic analysis; the quoted 
[Fe/H] errors include the uncertainties of the photometric 
indices and the published spread of the calibration. 



8. Absolute dimensions 

Absolute dimensions for BK Peg are presented in Tabic [TU1 
as calculated from the photometric and spectroscopic ele- 
ments given in Tables [7] and [5J As seen, masses and radii 
precise to 0.4-0.5% and 0.4-1.2%, respectively, have been 
established for the binary components. 

The V magnitudes and uvby indices for the compo- 
nents, included in Table [TUJ were calculated from the com- 
bined magnitude and indices of the system outside eclipses 
(Table [IJ and the luminosity ratios between the compo- 
nents (Table [7]). Within errors, the V magnitude and the 
uvby indices of the primary component agree with those 
measured at central secondary eclipse where 98% (y) of the 
light from the secondary component eclipsed (cf. Table [1} . 
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Table 10. Astrophysical data for BKPcg. 



Primary Secondary 



Absolute dimensions: 






ivl j M|?j 


1 ATA _1- n 007 
1.414 ± U.UU / 


l.zo < ± U.UUo 


z? / z? 

H/Hq 


l.yoo ± U.UUo 


1.4 f 4 zt U.Ul f 


log 9 (cgs) 


o.yyz ± U.UU4 


a oni i n m n 
4.zUl ± U.U1U 


■a /l — 1\ 

v sin z (kms ) 


16.6 ± 0.2 


13.4 ± 0.2 


b /i — 1\ 

Dsync (kms J 


18.3 ± 0.1 


13.6 ± 0.2 


Vpsync (kms ) 


18.3 ±0.1 


13.6 ± 0.2 


Vper/ 1 (kmS _1 ) 


18.5 ±0.1 


13.7 ±0.2 


Photometric data: 






V a 


10.473 ± 0.009 


11.080 ± 0.014 


(b-vT 


0.363 ± 0.007 


0.360 ± 0.008 


7721° 


1 4-4 + 01 3 


1 3Q + 01 6 


Cl° 


0.466 ± 0.014 


0.436 ± 0.017 


E(b-y) 


0.044 


±0.015 


T a 


6265 ± 85 


6320 ± 90 


M bo i 


2.90 ±0.06 


3.51 ±0.07 


log L/L 


0.74 ±0.02 


0.49 ±0.03 


B.C. 


-0.01 


-0.01 


M v 


2.91 ±0.06 


3.52 ±0.07 


Va - M v 


7.37 ±0.09 


7.37 ±0.10 


Distance (pc) 


298 ± 12 


298 ± 13 


Abundance: 






[Fe/H] 


-0.12 ±0.07 



a Observed rotational velocity 

b Equatorial velocity for synchronous rotation 

c Equatorial velocity for pseudo-synchronous rotation 

d Refers to periastron velocity 

e Not corrected for interstellar absorption/reddening 

Note: Bolometric corrections (BC) by Flower (1996 ) have been 
assumed, together with T e jjQ — 5780 K, BCq = —0.08, and 
M bolQ = 4.74. 



Table 11. Model parameters and average ages for BKPeg 
inferred from the observed masses and radii. 



Grid 


[Fc/H] 


Y 


Z 


Age (Gyr) 


Yonsei-Yale 


-0.190 


0.2540 


0.0120 


2.55 




-0.120 


0.2580 


0.0140 


2.70 




-0.105 


0.2590 


0.0145 


2.75 




-0.050 


0.2626 


0.0163 


2.85 


Victoria- Regina 


-0.190 


0.2629 


0.0125 


2.35 




-0.105 


0.2684 


0.0150 


2.50 




-0.039 


0.2735 


0.0173 


2.60 



suits, 6325 and 6345 K, within errors (Sect[7]). As seen from 
Table [S] (combined light), temperatures from the (b — y), c\ 
calibration by Alonso et al. (1996) and the b — y calibration 
by Ramirez & Melendez p005[) are in perfect agreement 
with that from the Holmberg et al. calibration. 

The measured rotational velocity for the secondary com- 
ponent is in perfect agreement with (pseudo) synchronous 
rotation, whereas the primary component seems to rotate 
at a slightly lower rate. 

The distance to BK Peg was calculated from the " clas- 
sical" relation (see e.g. CTB08), adopting the solar val- 
ues and bolometric corrections given in Table 1101 and 
A v /E(b - y) = 4.28 (Crawford & Mandwewala 1976). As 
seen, identical values are obtained for the two components, 
and the distance has been established to 4%, accounting 
for all error sources. Nearly the same distance is obtained 
if other BC scales, e.g. Code et al. (|1976p . Bessell et al. 
([TT)9"g|) and Girardi et al. ([21702]) . are adopted. Also, the em- 
pirical K surface brightness - T c g relation by Kervella et al. 
(|2004p leads to an identical and perhaps even more precise 
(2%) distance. We refer to Clausen (j2004|) and Southworth 
et al. (|2005p for details on the use of eclipsing binaries as 
standard candles. 



The interstellar reddening E(b - y) = 0.044 ± 0.015, 
also given in Table [T0l was determined from the calibration 
by Olsen (|1988p , using the uvbyf3 standard photometry for 
the combined light outside eclipses. Within errors, the same 
reddening is obtained from the indices observed during the 
central part of secondary eclipse. The new intrinsic-colour 
calibration by Karatas & Schuster (|2009|) leads to E(b — 
y) = 0.034. The model by Hakkila et al. (fi"997f yields a 
higher reddening of E(B - V) = 0.14 or E(b - y) = 0.10 
in the direction of and at the distance of BK Peg, whereas 
the maps by Burstein & Heiles f|1982[> and Schlegel et al. 
PW8"1) give total E(B - V) reddenings of 0.07 and 0.05, 
respectively. 

From the individual indices and the calibration by 
Holmberg et al. (2007J), we derive effective temperatures 
of 6265 ± 85 and 6285 ± 90 K for the primary and sec- 
ondary component, respectively, assuming the final [Fe/H] 
abundance. The temperature uncertainties include those of 
the uvby indices, E(b — y), [Fe/H], and the calibration itself. 
Compared to this, the empirical flux scale by Popper (1980) 
and the y flux ratio between the components (Tabled yield 
a well established temperature difference between the com- 
ponents of 55 ± 5 K (excluding possible errors of the scale 
itself). Consequently, we assign temperatures of 6265 and 
6320 K, which agree with the corrected spectroscopic re- 



9. Comparison with stellar models 

In the following, we compare the absolute dimensions ob- 
tained for BKPeg with properties of recent theoretical 
stellar evolutionary models. We have concentrated on the 
Yonsei-Yale (Y 2 ) grids by Demarque et al. ([2004 13 and 
the VRSS (scaled-solar abundances of the heavy elements) 
Victoria- Regina grids (VandenBerg et al., [20060 listed in 
Table ITTl The abundance, mass, and age interpolation rou- 
tines provided by the Y 2 group and the isochrone interpo- 
lation routines provided by the Victoria- Regina group have 
been applied. A summary of the Y 2 and VRSS grids and 
their input physics is given by CTB08. Here we just recall 
the following: The Y 2 models include He diffusion, whereas 
diffusion processes are not included in the VRSS models. 
The Y 2 models adopt the enrichment law Y = 0.23 + 2Z 
together with the solar mixture by Grevesse et al. (1996), 
and the VRSS models Y = 23544 + 2.2Z and the solar 
mixture by Grevesse & Noels ([1993). 

Both grids include core overshoot. For the Y 2 models, 
Aos = A ov /-ffp depends on mass and also takes into ac- 



http : //www . astro . yale . edu/demarque/yystar . html 
7 http : //wwwl . cadc-ccda.hia-iha.nr c-cnrc .gc . ca/cvo/| 
community/Viet or iaReginaModels/ 
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Fig. 5. Comparison between Y 2 (black) and VRSS (red) 
models for [Fe/H] = -0.105. Tracks (solid lines) for 1.2, 1.3, 
and 1.4 Mq and isochrones (dashed lines) for 2.0, 2.5, and 
3.0 Gyr are included. The 1.26 and 1.41 Mq components 
of BK Peg are shown. 




1.25 1.30 1.35 1.40 1.45 

M/M 



Fig. 6. Comparison between Y 2 (solid lines, black) and 
VRSS (dashed lines, red) models for [Fe/H] = -0.105. 
Isochrones for 1.0, 2.0, 2.5 (VRSS only), 2.75 (Y 2 only), 
3.0, and 4.0 Gyr are included. The 1.26 and 1.41 Mq com- 
ponents of BK Peg are shown. 



count the composition dependence of M™ ; " v 0. The ramp- 
ing of A os = X ov /H p from at M c c °£ v to 0.2 at 1.2 x M c c °£ v 
and higher masses is done in steps of 0.05. For the metallic- 
ity of BK Peg, M™£ v is 1 .2 Mq . The Victoria-Regina group 
adopts a somewhat different physically based description of 
the core overshoot with the free parameter F ovor depending 
on mass and metallicity and calibrated obscrvationally via 



Defined as "the mass above which stars continue to have 
a substantial convective core even after the end of the pre-MS 
phase." 




7000 6800 6600 6400 6200 6000 
T e „ [K] 

Fig. 7. BKPeg compared with Y 2 models for [Fe/H] = 
—0.12. Tracks for the component masses (solid lines) and 
isochrones for 0.5 and 1.0-4.0 Gyr (dashed; step 1.0 Gyr) 
are shown. Tracks (dotted) for [Fe/H] = —0.05 are included. 



cluster CMDs. For the metallicity of BK Peg, core convec- 
tion sets in around 1.15 Mq, and F over gradually increases 
to 0.55 at 1.70 Mq and is kept at this value for higher 
masses. Thus, for both grids, the 1.26 and 1.41 Mq com- 
ponents of BK Peg are in their ramping regions. 

In Fig. [SJ we compare Y 2 and VRSS main-sequence evo- 
lutionary tracks and isochrones for [Fe/H] = —0.1053 in the 
mass region of BKPeg. For 1.3 Mq the tracks nearly co- 
incide, whereas the TAMS positions differ for both lower 
and higher masses, probably related to the individual core 
overshoot recipes. In addition, somewhat different ages are 
predicted along the main-sequence tracks. This is also illus- 
trated by the isochrones in the mass-radius diagram shown 
in Fig. [51 which represents the most direct comparison with 
BK Peg, since the observed masses and radii are scale inde- 
pendent. As seen, similar ages of 2.75 (Y 2 ) and 2.50 (VRSS) 
Gyr are predicted for the components, with a slight pref- 
erence for the shape of the VRSS isochrone. Within the 
observed metal abundance range, [Fe/H] = —0.12 ± 0.07, 
the best fitting isochrones reproduce BK Peg nearly equally 
well; the corresponding average ages are listed in Table [Til 

Y 2 tracks for the observed masses and abundance of 
BK Peg are shown in Fig.[7J They are, like the VRSS tracks, 
slightly hotter. We also recall that the temperatures derived 
as part of the spectroscopic abundance analysis are slightly 
higher than those adopted from the uvby indices fTable [TU|) . 

In conclusion, the Y 2 and Victoria-Regina VRSS mod- 
els, including their core overshoot treatment and ramping 
recipes, are able to reproduce the observed properties of 
BK Peg. We note, however, that both grids predict a slightly 
lower age for the 1.41 Mq primary component than for the 
1.26 Mq secondary In the following section, we will check 
whether this tendency is seen for other well-studied binaries 
with similar component masses. 



9 VRSS models for the observed [Fe/H] = —0.12 are not avail- 
able 
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Table 12. Masses, radii, and abundances from Torres et al. ((2D09) for a subset of well-studied binaries with both 
components in the 1.15-1.70 Mq interval; see the text for details. 



System 


Colour 

r lgs. \v] I1UI 


M 

IVl pyx 

( M Q) 


IT) \ 


M sec 

fM \ 
\ M Q) 


R, 


fFc/Hl 
observed 


fFe/Hl 

V 2 
1 


Aee ( Gvrl 
V 2 


fFe /HI 

VP 


Aee ( Gvr) 

VP QQ 


GX Gem 


cyan 


1 488 


2 326 


1 467 


2 236 




0.00 


2.79 


0.000 


2.55 




±11 


±12 


±10 


±12 












BKPeg 




1 414 


1 987 


1 257 


1 473 


—0 12 


—0.12 


2.70 


—0.105 


2.50 






' ±7 


' ±8 


' ±5 


±17 


±7 










RW Aor 


TYlfl (XPTlt fl 
1 1 1 ci ti i... 1 1 L 1 a 


1 377 


1 786 


1 479 


2 062 


—0 07 


—0.08 


2.43 


—0.105 


2.15 




±21 


±43 


±19 


±44 


±11 










V442 Cvs 


red 


1 560 


2 073 


1 407 


1 663 




—0.10 


1.77 


—0.105 


1.65 




±24 


±34 


±23 


±33 












AD Boo 


green 


1.414 


1.612 


1.209 


1.216 


±0.10 


±0.10 


1.77 


±0.136 


1.50 






±9 


±14 


±6 


±10 


±15 










FS Mon 


orange 


1.632 


2.052 


1.462 


1.629 




±0.23 


1.44 


±0.226 


1.23 




±10 


±12 


±10 


±10 












VZHya 


brown 


1.271 


1.314 


1.146 


1.113 


-0.20 


-0.20 


1.20 


-0.190 


0.90 




±9 


±5 


±6 


±7 


±12 










V570 Per 


blue 


1.447 


1.521 


1.347 


1.386 


±0.02 


±0.02 


1.00 


0.000 


0.75 






±9 


±34 


±8 


±19 


±3 










HD71636 


gray 


1.514 


1.571 


1.288 


1.363 




0.00 


0.80/1.70 


0.000 


0.70/1.30 




±7 


±26 


±6 


±26 









10. Comparison with other binaries 

The review on accurate masses and radii of normal stars 
by Torres et al. (|2009|) lists 20 binarie£3 with both compo- 
nents in the 1.15-1.70 Mq interval where core overshoot is 
ramped up in the Y 2 and Victoria-Regina models. Here, we 
will concentrate on the subset listed in TableHU where most 
of the nearly equal-mass binaries have been excluded. Eight 
binaries have been selected because the masses of their com- 
ponents differ significantly; we have adopted a lower mass 
difference limit of 0.10 Mq. In addition, we have included 
GXGem, since its components have evolved close to the 
TAMS region. As mentioned in Sect. [Q Lacy et al. ([2"UUg|) 
found that for GX Gem, the lowest a ov consistent with ob- 
servations is approximately 0.18. [Fe/H] has been measured 
for five of the systems; we refer to Appendix [A~l for details 
on the abundance analysis of BW Aqr. 

Fig. [5] shows the nine binaries together with Y 2 evo- 
lutionary tracks for the measured masses. For GXGem, 
V442Cyg, FS Mon, and HD 71636, [Fe/H] has not been 
measured. The values listed in Table [T^] were adopted only 
because they lead to reasonable model fits for the ob- 
served effective temperatures. The isochrones correspond 
to the average age inferred from the measured masses and 
radii. To illustrate the abundance dependence, evolution- 
ary tracks and isochrones for solar metallicity have been 
included for V442 Cyg and FS Mon. As seen, the nine bina- 
ries cover the main-sequence almost from the ZAMS to the 
TAMS. 

In Figs. [9] (Y 2 ) and [TUl ( VRSS). we concentrate on com- 
parisons based on the scale-independent masses and radii. 
Since abundance interpolation software is not available for 
the VRSS models, we have used those with [Fe/H] clos- 
est to the measured values. The average ages are listed in 
Table [T^l and we note the following: 

— The VRSS models predict lower ages than the Y 2 mod- 
els for all systems. 



The 20 binaries can now be supplemented by BKPeg and 
V1130 Tau (Clausen et al. |20T0]> . 



— For the four binaries with components in the upper half 
of the main-sequence band (GX Gem, BK Peg, BW Aqr, 
V442 Cyg) , both grids tend to predict lower ages for the 
more massive component than for the less massive one, 
although the difference is certainly very marginal for the 
nearly equal-mass system GX Gem. These systems span 
the 1.26-1.56 Mq interval. 

— The Y 2 isochrones fit FSMon (1.63 ± 1.46 Mq) bet- 
ter than the VRSS isochrones, which predict a higher 
age for the more massive component than for the less 
massive one. 

— Both grids fit the less evolved binaries AD Boo and 
VZ Hya quite well, perhaps with a slight preference for 
Y 2 . 

— For the little evolved binary V570Per, VRSS is 
marginally better than Y 2 . 

— HD 71636 is not fitted well by any of the grids. The very 
different ages predicted for the components are listed 
in Table [T2] We believe, however, that this is due to 
problems with the published radii. A re-analysis of the 
light curves by Henry et al. ([2006) indicates that the 
ratio between the (relative) radii is sig nificantly more 
uncertain than given by the authorfH. most probably 
because the secondary eclipse is far from being well- 
covered. 

It is outside the scope of this paper to propose specific 
modifications of the model physics, but we believe that this 
sample of binaries (except HD 71636) can be used to fine- 
tune the core overshoot treatment in terms of mass and 
metal abundance. It might also be relevant to include in- 
vestigations of other model ingredients, e.g. diffusion pro- 
cesses and the adopted helium-to-metal enrichment ratio, 
and to clarify why the Y 2 models predict higher ages than 
the VRSS models. In order to avoid possible interpolation 
errors, specific models for the observed masses and metal 
abundances should be calculated and small age steps ap- 
plied. On the observational side, spectroscopic metal abun- 

We have used the larger errors listed by Torres et al. (2009), 
which also seem to be underestimated. 
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Fig. 8. Y 2 evolutionary tracks (full drawn lines) for the binaries included in Table Q21 The isochrones (dashed lines) 
correspond to the average ages inferred from the masses and radii for the [Fe/H] listed in the table. For HD 71636, the 0.80 
Gyr isochrone is plotted. For V442Cyg and FSMon, the dotted tracks and isochrones (1.83 and 1.38 Gyr, respectively) 
represent solar abundance models. Component masses (in units of Mq) are given in the lower left corners. 



dances should be determined for the four binaries lacking 
this information, and it should perhaps be considered to 
re-observe BWAqr, V442 Cyg, and HD 71636 to improve 
their masses and/or radii. 

Also, we suggest to supplement the sample by the 
unique K0IV+F7V binary AI Phe (1.23 + 1.19 Mq), re- 
cently re-discussed by Torres et al. (I2009p . as well as by a 
number of new F-type systems we are presently studying. 
In addition, it would be highly relevant to include F-type 
binary members of open clusters. 

11. Summary and conclusions 

From state-of-the-art observations and analyses, precise 
(0.4-1.2%) absolute dimensions have been established for 
the components of the late F-type detached eclipsing bi- 
nary BKPeg (P = 5?49, e = 0.0053); see Table HI A 
detailed spectroscopic analysis yields an iron abundance 
relative to the Sun of [Fe/H] = -0.12 ± 0.07 and similar 
relative abundances for Si, Ca, Sc, Ti, Cr, and Ni. The 
measured rotational velocities are 16.6 ±0.2 (primary) and 
13.4±0.2 (secondary) kms -1 . For the secondary component 
this corresponds to (pseudo)synchronous rotation, whereas 
the primary component seems to rotate at a slightly lower 
rate. 

The 1.41 and 1.26 Mq components of BKPeg have 
evolved to the upper half of the main-sequence band. 
Yonsei-Yale and Victoria-Regina solar scaled evolution- 



ary models for the observed metal abundance reproduce 
BKPeg at ages of 2.75 and 2.50 Gyr, respectively, but tend 
to predict a lower age for the more massive primary com- 
ponent than for the secondary. If real, this might be due 
to less than perfect calibration of the amount of convective 
core overshoot of the models as function of mass (and metal 
abundance) . 

For this reason, we have performed model comparisons 
for a sample of eight additional well-studied binaries with 
component masses in the 1.15-1.70 Mq interval where con- 
vective core overshoot is gradually ramped up in the mod- 
els; see Table [T2J We find that a) the Yonsei-Yale models 
systematically predict higher ages than the Victoria-Regina 
models, and that b) the three other most evolved systems in 
the sample share the age difference trend seen for BK Peg. 

We propose to use the sample to fine-tune the core over- 
shoot treatment, as well as other model ingredients, and to 
clarify why the two model grids predict different ages. The 
sample should be expanded by a number of new F-type sys- 
tems under study, binary cluster members, and the unique 
K0IV+F7V binary AlPhe (1.23+1.19 Mq). 
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Fig. 9. Comparison between Y 2 isochrones (full drawn 
lines) and the binaries listed in Table H21 we refer to the 
table for colour codes, [Fe/H] and ages. For V442 Cyg and 
FS Mon, the dashed lines are solar abundance isochrones 
for 1.83 and 1.38 Gyr, respectively. For HD 71636, the 0.80 
Gyr isochrone is plotted. 




Fig. 10. Comparison between VRSS isochrones (full drawn 
lines) and the binaries listed in Table [T^l we refer to the 
table for colour codes, [Fe/H] and ages. For V442 Cyg and 
FS Mon, the dashed lines are solar abundance isochrones 
for 1.67 and 1.27 Gyr, respectively. For HD 71636, the 0.70 
Gyr isochrone is plotted. 
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Appendix A: Chemical abundances for BWAqr 

Absolute dimensions for the late F-type eclipsing binary 
BWAqr were published by Clausen (|199ip . who mentioned 
that the uvby/3 photometry indicates a metallicity slightly 
above solar. We note that the primary component (star A 
in Clausen 11 99 Ij) , eclipsed at the deeper minimum at phase 
0.0, is less massive, smaller, but hotter than the secondary 
component. 

Here we present the results from an abundance anal- 
ysis based on two high-resolution spectra observed with 
the FEROS fibre echelle spectrograph at ESO, La Silla in 
August 1999; see Table IA.1I Details on the spectrograph, 
the reduction of the spectra, and the basic approach fol- 
lowed in the abundance analysis are described by CTB08. 
As for BKPeg, we have used the VWA tool for the abun- 
dance analysis, and we refer to Sect. [7] for further infor- 
mation on atmosphere models, atomic data information, 
log(<?/) adjustments etc. One important difference is, how- 
ever, that for BW Aqr disentangling is not possible, and the 
analysis is therefore based on double-lined spectra. 



Table A.l. Log of the FEROS observations of BWAqr. 



hjd-2 4oo oorr 


phase 


t 6 


S/N^ 


51393.82092 


0.1290 


3600 


140 


51394.70920 


0.2612 


3600 


200 



" Refers to mid-exposure 
b Exposure time in seconds 

c Signal-to-noise ratio measured around 6070 A 

Table A. 2. Astrophysical data adopted for the abundance 
analysis of BWAqr. 



Primary Secondary 



^eff (K) 


6300 


6225 


log g (cgs) 


4.075 


3.981 


vs'm i (km s _1 ) 


14.0 


14.0 


Vmicro (km S ) 


1.70 


1.50 



The effective temperatures, surface gravities, rota- 
tional velocities, and microturbulence velocities listed in 
Table IA.2I were adopted. The temperatures were deter- 
mined by requiring that Fe I abundances were independent 
of line excitation potentials. They are slightly lower than 
determined by Clausen (|199ip . who obtained 6450 ± 100 
K (primary) and 6350 ± 100 K (secondary). However, new 
and better (unpublished) uvby photometry for BW Aqr and 
the calibration by Holmberg et al. (|2007p lead to 100 K 
lower values, assuming a reddening of E(b — y) — 0.03. 
Microturbulence velocities were tuned until Fe I abun- 
dances are independent of line equivalent widths. 

The abundances derived from all useful lines in both 
spectra are presented in Table IA.3I The equivalent widths 
measured in the two double-lined spectra are listed in 
Tables 16 (primary) and 17 (secondary), which will only 
be available in electronic form. Comparing the results from 
the two spectra, we find that they agree within 0.05 dex. We 
have only included lines with measured equivalent widths 
above 10 mA and below 45 mA (primary) and 55 mA 
(secondary). The lines are diluted by factors of about 2.2 
(primary) and 1.8 (secondary), meaning that lines with in- 
trinsic strengths above 100 mA arc excluded. The [Fe/H] 
results for the two components differ by 0.06 dex, but 
for each component the results from Fe I and Fe n lines 
agree well. The mean value for all measured Fe lines is 
[Fe/H] =-0.07 ± 0.01 (mis of mean). 

Changing the model temperatures by ±100 K modifies 
[Fe/H] from the Fe I lines by about ±0.10 dex, whereas 
almost no effect is seen for the Fe n lines. If 0.30 kms -1 
higher microturbulence velocities are adopted, [Fe/H] de- 
creases by about 0.05 dex for both neutral and ionized lines. 
Taking these contributions to the uncertainties into ac- 
count, we adopt [Fe/H] = -0.07±0.11 for BWAqr. Similar 
abundances are obtained for the few other ions listed in 
Table IA~31 

List of Objects 

'BWAqr' on page U\ 
'GX Gem' on pageQ] 
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Table A. 3. Abundances ([El./H]) for the primary and 
secondary components of B W Aqr determined from the two 
FEROS spectra. 



Primary Secondary 



Ion 


[El./H] 


rms 


N t /N; 


[El./H] 


rms 


N t /N, 


Si I 


-0.12 


0.02 


3/2 


-0.05 


0.12 


13/7 


Cal 


-0.18 


0.13 


3/2 


-0.00 


0.09 


7/4 


Sell 








-0.08 


0.08 


3/2 


Til 


-0.02 


0.09 


3/2 








Cri 


-0.12 


0.17 


4/2 


-0.04 


0.12 


5/4 


Fel 


-0.11 


0.16 


65/41 


-0.05 


0.11 


80/50 


Fen 


-0.09 


0.18 


14/8 


-0.03 


0.11 


14/8 


Nil 


-0.05 


0.14 


18/11 


-0.12 


0.10 


15/10 



Note: Nj is the total number of lines used per ion, and N; is 
the number of different lines used per ion. Ions with at least 3 
lines measured are included. 
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